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information encoded in the dispersion relation for the D-function, the effects due 
to the interrelation between spacelike and timclike domains, and the effects due to 
the nonvanishing pion mass. The latter plays a crucial role in this analysis, forcing 
the Adler function to vanish in the infrared limit. Within the developed approach 
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additional input. The obtained result is found to be in reasonable agreement with 
the experimental prediction for the Adler function in the entire range of momenta 
< Q 2 < co. 
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1. Introduction 

The hadronic vacuum polarization function II(g 2 ) and the Adler function D(Q 2 ) are 
physical observables of central importance in various fields of elementary particle 
physics. In particular, these quantities are essential for the analysis of strong 
interaction processes such as electron-positron annihilation into hadrons PJ El El H] 
and r lepton decay [3 EJ- Furthermore, the Adler function plays an important role 
when confronting precise experimental measurements of some electroweak observables 
with their theoretical predictions, giving rise to decisive tests of the Standard Model 
and furnishing stringent constraints on possible new physics beyond it. Specifically, due 
to the hadronic excitations in the photon vacuum polarization, the muon anomalous 
magnetic moment [HE] an d the shift of the electromagnetic fine structure constant P1ITU] 
receive appreciable contributions from the strong interactions. The latter constitute in 
fact the main source of theoretical uncertainty when computing these observables, since 
the physical energy scales involved are such that the nonperturbative effects begin to 
be relevant. 

To date, there is no systematic method for calculating the Adler function in the 
whole range of energies < Q 2 < oo. Nonetheless, in the asymptotic ultraviolet region 
D(Q 2 ) may be approximated by the power series in the strong running coupling a s (Q 2 ) 
by making use of perturbation theory. However, unphysical singularities of a s (Q 2 ) 
(e.g., the one-loop Landau pole), being artifacts of the perturbative computations, 
invalidate this approach in the infrared domain. In turn, this significantly complicates 
the theoretical description of low-energy experimental data, and eventually forces one 
to resort to various models and phenomenologically inspired approximations. 

The "semi-experimental" method for obtaining the Adler function in the infrared 
domain is to actually employ the data on the i?-ratio of e + e~ annihilation into hadrons, 
in conjunction with the dispersion relation [2]. Specifically, one first constructs R(s) by 
merging its low-energy experimental data with its high-energy perturbative prediction, 
and then integrates the dispersion relation for the Adler function |2J. This way of 
modeling of the i?-ratio entails certain complications. In particular, one has to properly 
take into account the effects due to the analytical continuation of the perturbative 
results for D(Q 2 ) into the timelike domain (such as the resummation of the so-called 
7r 2 -terms, see Ref. JI]), and resolve the ambiguities related to the matching condition 
between the experimental and perturbative inputs for R(s). In addition, the e + e~ data 
are not always suitable for the extraction of the Adler function, basically due to large 
systematic uncertainties in the infrared domain. Nevertheless, recent progress in the 
study of the r lepton decays offers a way to overcome the latter difficulty Specifically, 
for the energies below the mass of the r lepton, the e + e~ data can be substituted 
(up to isospin breaking effects [12]) by a precise inclusive vector spectral function [13j, 
extracted from the hadronic r decays. It is worthwhile to note also that some insights 
on the infrared behavior of D(Q 2 ) may be gained from chiral perturbation theory (see, 
e.g., paper [JJ] and references therein) and lattice simulations [To] . 



3 



An important source of the nonperturbative information about the hadron 
dynamics at low energies is provided by the relevant dispersion relations. The latter, 
being based on the general principles of the local Quantum Field Theory (QFT), 
supply one with the definite analytic properties in a kinematic variable of a physical 
quantity at hand. The idea of employing this information together with the perturbative 
treatment of the renormalization group (RG) method forms the underlying concept of 
the so-called "analytic approach" to QFT. It was first proposed in the framework of 
the Quantum Electrodynamics and applied to the study of the invariant charge of 
the theory [16 a . Later on, it was argued (see, e.g., papers |U [T71 UHl HHJ EH! and 
references therein) that a similar method can also be useful for studying the non- 
Abelian theories. Eventually, proceeding from these motivations, the analytic approach 
to Quantum Chromodynamics (QCD) has been developed j2Qj. Some of the main 
advantages of this method are the absence of the unphysical singularities and a fairly 
good higher-loop and scheme stability of outcoming results. The analytic approach 
has been successfully employed in studies of the strong running coupling [2H1 EI], 
perturbative series for QCD observables (the so-called "Analytic perturbation theory" 
(APT), see papers |221 ESI EH EE] and references therein), hadron spectrum [2E], pion 
form factor and some other intrinsically nonperturbative aspects of the strong 
interaction [2D EE]- 

The primary objective of this paper is to derive novel integral representations for 
the Adler function and the i?-ratio in a general framework of the analytic approach to 
QCD. These representations contain the nonperturbative information captured in the 
dispersion relation for D{Q 2 ), the effects due to the interrelation between spacelike and 
timelike domains, and the effects due to the nonvanishing pion mass. In addition, we 
compute the Adler function within the approach at hand, by employing its perturbative 
approximation as the only additional input, and compare the obtained result with the 
relevant experimental data. 

The layout of the paper is as follows. In Sec. El the dispersion relation for the 
Adler function D(Q 2 ) and its relation to the .R-ratio of the e + e~ annihilation into 
hadrons are briefly reviewed. In Sec.|3]new integral representations for D(Q 2 ) and R(s), 
which involve a common spectral function, are derived. The calculation of the Adler 
function within the developed approach, its comparison with the experimental prediction 
for D(Q 2 ), and a discussion of the obtained results are presented in Sec. |U 

2. Dispersion relation for the Adler function 

Let us consider the process of the electron-positron annihilation into hadrons in the 
lowest order of the electromagnetic interaction (see Fig.HJ. The corresponding Feynman 
amplitude is given by 

v(p 2 )eru( Pl )\(T\J,(q)\0), (1) 
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Figure 1. The Feynman diagram for the electron-positron annihilation into hadrons. 

where u(pi) and ^(^2) are the Dirac spinors of the electron and positron respectively, 
s = q 2 = (pi + P2) 2 is the timelike momentum transferred, T denotes a final hadron 
state, and J^q) stands for the quark current. The probability of the transition of an 
electron-positron pair into hadrons is proportional to the square of the amplitude (JTJ) 
summed over the states T. The hadronic part of the latter can be represented as 

^M) = E <° M-dM (r | J v { q )\ o> = (o I J M (-g) Mq)\ o> , (2) 

r 

where the completeness ^rl^X-H = 1 nas been employed. It is worth emphasizing 
here that A M ^(g 2 ) exists only for q 2 > 4m 2 , since otherwise no hadron state T could be 
excited. In particular, A^ u (q 2 ) vanishes identically below the two-pion threshold (i.e., 
for q 2 < 4m 2 ), see Ref. |T] for the details. 

The hadronic tensor A /iy (g 2 ) is proportional to the discontinuity of the function 

IV(g 2 ) = i j Ae"' 1 (0 |T{J M (x) J„(0)}| 0) = (q^ - g^q 2 )Tl(q 2 ) (3) 

across the physical cut q 2 > 4m 2 along the positive semiaxis of real q 2 . In Eq. (|3*|) 
Il(g 2 ) stands for the hadronic vacuum polarization function. The latter satisfies the 
once-subtracted dispersion relation QjJ |2| 

ntf) = m^-( g »-^ j T (t _gg_ y) rf ., (4) 

where m w = (134.9766 ± 0.0006) MeV [22] is the mass of the vr meson. In Eq. (gj) R(s) 
denotes the measurable ratio of two cross-sections 

7-,/ x It n-r/ • \ tt/ ■ m °" (e + e~ — > hadrons; s) 

R(s) = — hm n(s - 2£ - n (s + ze = -i— ^, 5 

with s being the center-of-mass energy squared of the annihilation process. It is worth 
noting here that R(s) = for s < 4m 2 , because of the kinematic restrictions mentioned 
above. 

For practical purposes it proves convenient to deal with the Adler function j2] , which 
is defined as the logarithmic derivative of the hadronic vacuum polarization function (JHJ 
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and, therefore, does not depend on the choice of subtraction point s' in dispersion 
relation Pjl. In Eq. © Q 2 = —q 2 > denotes a spacelike momentum. The Adler 
function (jHJ plays an indispensable role for the simultaneous processing of the timelike 
and spacelike experimental data, see Refs. [2 El • Indeed, since the perturbation theory 
and the RG method are not applicable directly to the study of the observables depending 
on the timelike kinematic variable, for the self-consistent description of the latter one 
first has to relate the timelike experimental data with the perturbative results. Here, 
the required link between the i?-ratio (jSJ) and the Adler function (jHJ) can be obtained 
by differentiating Eq. that results in the dispersion relation |5j 



As it has been mentioned in the Introduction, this equation is commonly employed 
for restoring D(Q 2 ) from the e + e~ experimental data. On the other hand, the inverse 
form of relation ([7]). which enables one to continue an explicit expression for the Adler 
function into timelike domain, can be derived by integrating Eq. (jUJ) between finite 
limits. Ultimately, this leads to 



where the integration contour lies in the region of analyticity of the integrand, see 
Ref. jllj . In particular, Eq. (jSJ) provides the only known way for obtaining an explicit 
expression for R(s), using theoretical prediction for D{Q 2 ) as input (see also Refs. 
for details). 

3. Novel integral representations for D(Q 2 ) and R(s) 

Before proceeding to the derivation of the integral representations for D(Q 2 ) and R(s), 
it is worth emphasizing that the dispersion relation (J7|) imposes stringent constraints 
on the form of the Adler function. Specifically, Eq. (JZj) implies that D(Q 2 ) is an 
analytic function in the complex Q 2 -plane with the only cut Q 2 < —Am 2 along 
the negative semiaxis of the real Q 2 . In addition, given that (i) R(s), being a 
physical quantity, assumes finite values, and (ii) its asymptotic ultraviolet behavior 
R(s) ~ l + (9[ln~ 1 (s/A 2 )] when s — > oo, one concludes from Eq. (J7J) that D(Q 2 ) vanishes 
in the infrared limit Q 2 —>■ 0. 

Let us for the moment turn off the strong interactions, and neglect all effects due to 
quark masses (these latter effects will be disregarded throughout the paper, even when 
the strong interaction are turned back on; a brief discussion of their possible relevance 
will be given in the next section). In this case, the i?-ratio © is determined by the 
parton model prediction PQ: 



where the overall factor N c Q 2 is omitted throughout, N c = 3 is the number of 
colors, Qf denotes the charge of the quark of the fth flavor, and 8(x) stands for the 




(7) 




(8) 



R {s) = 6{s - Ami) 



s > 0, 



(9) 
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Figure 2. The integration contour in Eq. (JHJ) for s > 4mJ. The physical cut of the 
Adler function D(—Q is shown along the positive semiaxis of real £. 



Heaviside step-function. In turn, as it follows from the dispersion relation (J7J), Eq. Q 
corresponds to the following zeroth order prediction for the Adler function: 

d °^ = q^< q2>0 - (10) 

Obviously, this expression vanishes in the infrared limit Q 2 — > 0, and satisfies the 
physical condition Dq(Q 2 ) — > 1 when Q 2 — > oo. Thus, along with the "standard" 
intrinsically nonperturbative contributions j^Ol EI] the Adler function also receives 
power corrections due to the nonvanishing pion mass, which turn out to be important 
for Q < 2 GeV: 

^i^+EL^)"- ^ 2 > 4 ^- (id 

It is worth noting that the modification of the parton model prediction (jl(Jj) due to 
the mass of the 7r meson is crucial for arriving at a form of D(Q 2 ) which agrees 
with the low-energy experimental data. Had one neglected the pion mass, one would 
have instead obtained as a zeroth order approximation Dq(Q 2 ) = 1, which would have 
been inconsistent with the infrared behavior of the Adler function extracted from the 
experiment (see Fig. EJ). Remarkable as it may seem, to the best of our knowledge, 
Eq. (|1U|) does not appear anywhere in the existing literature. 

The next step is to turn the strong interactions back on, and derive new integral 
representations for R(s) and D(Q 2 ), which involve a common spectral function, to be 
denoted by p D (o")- The inclusion of the effects due to the strong interaction modify 
Eqs. P) and (fT0|) : 

R(s) =0(s-4ml)+r(s), (12) 
D « 2 > = QOT + rfW2) - (13) 
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As it has been mentioned above, the ultraviolet behavior of the strong correction d(Q 2 ) 
can be computed in the framework of perturbation theory, and the corresponding r(s) 
may be obtained from it through Eq. (jSJ. Assuming that an explicit "exact" expression 
for the Adler function (jl3j) is available, one can restore R(s) (|12|) by making use of 
relation (|8|). For the energies below the two-pion threshold (i.e., for < s < Am 2 ) the 
integration of Eq. (JSJ) leads to R(s) = 0, in conformity with the kinematic restrictions 
mentioned above. For the energies above the two-pion threshold it is convenient to 
choose the integration contour in Eq. (jHJ) in the form presented in Fig. |2l since the strong 
correction d(Q 2 ) vanishes at the ultraviolet asymptotic. In this case the only nontrivial 
contribution into the integral along the circle of the infinitely large radius comes from 
the parton model prediction (jlUj) . whereas the strong correction d(Q 2 ) contributes only 
to the integral along the cut of D(Q 2 ). Eventually, this results in 



R(s) = 6(s - Am 



2> 





r ( \ da 


1 + 


/ Pd(<t) — 







(14) 



where the spectral function is determined as the discontinuity of the Adler function 
across the cut 

/°d(c) = 7T- hm [D(-a - is) - D(-a + is)] , a > Am\. (15) 

At the same time, p D (a) (|T3|) can also be related to R(s), by differentiating Eq. (jSj) with 
respect to In s: 

PM = o > Ami (16) 

a ma 

that allows one to extract p D (c) from the relevant experimental data. In turn, the Adler 
function can be represented in terms of the spectral function p u (cr) by integrating the 
dispersion relation (J7J) with R(s) given by Eq. ()14|) . Carrying out the integration by 
parts, one arrives at 

[ f°° . S a-Aml da' 



D(Q 2 ) Q 



Q2 + 4m 2 



(17) 



Umi <y + Q 2 o 

It is worth noting that in deriving the integral representations (|T4*j) and (JT7j) we 
have employed only Eqs. (jZJ) and (jHJ), the parton model prediction and the fact 
that the strong correction d(Q 2 ) vanishes in the asymptotic ultraviolet limit Q 2 — > 
oo; no additional approximations nor model-dependent assumptions were involved. 
Similarly to the perturbative approach, the QCD scale parameter A remains the basic 
characterizing quantity of the theory. Besides, all the effects due to the interrelation 
between the spacelike and timelike observables are automatically accounted for by the 
representations ()14j) and (fTTj) . It is important to mention also that in the limit of the 
vanishing pion mass m^ = the derived relations (|T4*j) and (fTTj) become identical to 
those obtained within the massless APT [22J, namely 



Pd(*) — , (18) 

^ 1 r Pn{°) 



D APT (Q 2 ) = 1+ / J2 rh^da. (19) 
a + Q 2 
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4. Discussion and conclusions 

We hasten to emphasize that the new integral representation for the Adler function (fTTj) 
incorporates the same nonperturbative constraints on D(Q 2 ) as those contained in the 
dispersion relation ((ZJ); namely, D{Q 2 ) of (fTTj) possesses correct analytic properties 
in the Q 2 variable, and vanishes in the infrared limit Q 2 — > 0. However, D(Q 2 ) 
of (fTTj) has the advantage of being expressed as the discontinuity of itself across the 
cut (j!5)) . This fact brings about certain simplifications in the theoretical analysis of 
the Adler function. Specifically, as it has been noted in the Introduction, the standard 
extraction of D{Q 2 ) through Eq. constitutes a two-step procedure: first, one has to 
construct R{s) for the entire energy range 4m 2 < s < oo (encountering the complications 
mentioned above), and only then integrate Eq. ((7|). Instead, Eq. (|T7|) eliminates the 
intermediate step of building up R(s), allowing the reconstruction of D(Q 2 ) through 
the spectral function (jl5|) . At the same time, one is still able to incorporate the 
relevant experimental information about R{s) into the spectral function p v (a) by virtue 
of Eq. ()16|). that, however, turns out to be a bit more complicated technically, since the 
numerical differentiation of the experimental data on R(s) is required. 

In order to exploit Eq. (fTTj) one should furnish an input for the central quantity 
of the approach at hand, namely, the spectral function p D (cr) (fTS|) . In general, p D (a) is 
supposed to embody all available information about functions D(Q 2 ) and R(s). In this 
paper we restrict ourselves to the study of the perturbative contributions to p D ((x) only. 
Following the same motivations as those of APT [22], the latter can be obtained by 
identifying D(Q 2 ) in Eq. (|T3jl with its perturbative approximation, which at the £-loop 
level reads (see paper [32] and references therein) 

dSAQ 2 ) = i + Yt j=l d i H e) (Q 2 )Y > Q 2 - oo. (20) 

It is worth noting that the resulting perturbative spectral function p pert (cx) receives no 
contributions from spurious singularities of a s (Q 2 ), since the former is determined as 
the discontinuity of D peTt (Q 2 ) across the physical cut. Specifically, at the one-loop level 
(di = 1/vr) one has PpitO) = (4/A))[hi 2 (cr/A 2 ) + tt 2 ]" 1 , where /3 = 11 - 2n f /3, and m 
is the number of active quarks. 

It turns out that, even with the one-loop perturbative approximation of the spectral 
function, Ppe It (cr), representation (fTTjl is capable of providing an output for D(Q 2 ) which 
is compatible with its experimental prediction. Indeed, as can be seen in Fig. El for the 
entire energy range, < Q 2 < oo, the result obtained is in good agreement with the 
experimental behavior of D(Q 2 ) extracted from the ALEPH data on the inclusive vector 
spectral function [13], in the way described above. The value of the scale parameter 
A = 335 MeV has been estimated for the case of rif = 2 active flavors by fitting 
D{Q 2 ) (fTTj) to its experimental prediction by means of the least-squares method. It is 
worth mentioning that the Adler function (jl7)l is stable with respect to the higher-loop 
perturbative corrections. Indeed, calculations up to four loops have revealed that the 
relative difference between the £-loop and the (£ + l)-loop expressions for D{Q 2 ) (|T7j) 
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Figure 3. Comparison of the Adler function 117(1 corresponding to p^rti' 7 ) (solid 
curve) with its experimental behavior (•) extracted from |13| . The zeroth order 
prediction for D(Q 2 ) (|10|l is shown by the dashed curve, the massless APT case <|19|) 
is denoted by the dot-dashed curve, and the dot-dot-dashed curve corresponds to the 
one-loop perturbative approximation (|20|l of the Adler function. 



is less than 4.9%, 1.5%, and 0.3%, for £ = 1, £ = 2, and £ = 3, respectively, for 
< Q 2 < oo (the estimation [33] of the four-loop expansion coefficient cLt has been 
adopted here). 

In order to illustrate the significance of the mass of the tt meson within this 
approach, it is worth presenting the massless APT prediction of D(Q 2 ) (fl"9"|) computed 
by making use of p^rti 17 ) (dot-dashed curve in Fig. [3J. In this case, one arrives at 
a result, which is free of infrared unphysical singularities (to be contrasted with the 
perturbative dot-dot-dashed curve in Fig. [3J, but fails to describe the experimental 
behavior of the Adler function in the low-energy domain Q < 1 GeV, where the 
effects due to the pion mass become appreciable. Of course, in the framework of the 
massless APT, the infrared behavior of the Adler function can be greatly improved 
following the procedure introduced in [33]; this consists essentially in carrying out 
an appropriate resummation of threshold singularities, and introducing into ()18|) and 
(fT9"|) effects from nonperturbative light quark masses. The necessary nonperturbative 
information on the quark masses is furnished from the study of Schwinger-Dyson 
equations and quark condensates. Specifically, one extracts a particular momentum- 
dependence for the dynamical (effective) quark mass, which interpolates between the 
constituent and current mass values. We believe however that the method developed 
here presents certain advantages compared to that of [33] • This is so because, in addition 
to the possible ambiguities stemming from the truncation and general treatment of the 
Schwinger-Dyson equations, the final numerical implementation of the procedure of [33] 
seems to require rather elevated constituent masses for the light quarks {m u = ma = 
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250 MeV); instead, the only phenomenological parameter appearing in Eqs. (|14|) and (|17|) 
is the measurable mass of the pion. Finally, again in the framework of massless APT, 
the infrared behavior of the Adler function has been brought into qualitative agreement 
with experiment, but at the expense of resorting to additional assumptions, such as the 
vector meson dominance |H5| . 

As mentioned in the previous section, in this paper we have omitted effects due 
to the quark masses. Of course, this is not to say that such effects are not relevant. 
As a matter of fact, and in addition to their role in bringing the predictions of APT in 
agreement with experiment, as exemplified in [21], they could be of relevance even within 
the approach at hand. In particular, in the present analysis such effects can be retained 
by employing the perturbative prediction for the Adler function, which accounts for 
the quark masses (see paper jHHj and references therein), instead of D t (Q 2 ) given by 
Eq. (J20J) - We hope to present a detailed study of this issue in the near future. 

The results obtained in this paper can be further scrutinized through the 
incorporation of higher order perturbative contributions into the spectral function p v (a), 
along with the effects due to the quark masses, power corrections due to gluon and quark 
condensates [SOI EE] , and restrictions imposed by the low-energy experimental data. It 
would also be of apparent interest to apply the present approach to the estimation of 
the hadronic contributions to the muon anomalous magnetic moment and to the shift 
of the electromagnetic fine structure constant. 
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